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Abstract 


Given  the  traffic  flow  from  each  source  to 
each  destination  in  a network  and  given  the  aggre- 
gate traffic  in  each  link,  ji/e  want  ta  find  if  there 
is  any  looping  of  traffic;'  A careful  definition  of 
looping  shows  that  the  question  is  equivalent  to 
whether  some  of  the  aggregate  link  flows  cein  be  re- 
duced without  increasing  any  of  the  others.  It  is 
then  shown,  through  the  use  of  duality  in  linear 
programming,  that  an  aggregate  flow  is  loop  free 
^iff  all  the  traffic  follows  shortest  routes  for 
some  assignment  of  positive  lengths  to  the  links. 

It  is  further  shown  that  there  is  a finite  set 
of  these  length  assignments,  dependent  only  on  the 
topology  of  the  network,  such  that  every  shortest 
route  flow  is  a shortest  route  flow  for  one  of 
those  special  assignments.  Finally,  it  is  shown 
that  any  loopfree  flow  can  be  realized  by  a rout- 
ing in  which  the  sum,  over  all  destinations,  of  the 
number  of  alternate  route  links  required  to  reach 
that  destination,  is  at  most  the  number  of  links 
minus  the  number  of  nodes. 


Introduction 

When  the  data  in  a communication  netv;ork  oc- 
casionally travel  in  loops,  these  loops  generate 
both  an  unnecessary  loss  of  resources  and  also  an 
unnecessary  increase  in  delay.  Our  purpose  here  is 
to  define  a generalized  form  of  looping  and  to  pro- 
vide both  a mathematical  basis  and  heuristic  under- 
standing for  this  phenomena.  We  do  not  develop  any 
new  routing  algorithms  here,  but  the  results  pro- 
vide insight  into  the  weaknesses  of  current  algo- 
rithms and  into  potential  directions  for  overcom- 
ing these  weaknesses. 

Consider  a network  of  n nodes  and  L directed 
links.  We  denote  the  nodes  by  the  integers  1,..., 
n and  denote  the  links  by  integer  pairs,  (i,k)  de- 
noting a link  from  i to  k.  We  assume  that  the  net- 
work is  connected  in  the  sense  that  for  each  pair 
of  nodes  j _ there  is  some  directed  path  of  links. 
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say  ((i,k),(k,f),(a,m),...,  (h,j))  going  from  i to 
j.  Let  r^(j)  denote  the  time  average  data,  in  bits 

per  second,  entering  the  network  at  node  i and 
destined  for  some  receiver  at  node  j.  Similarly, 
let  ^ 0 denote  the  time  average,  in  bits 

per  second,  of  the  traffic  from  all  Inputs  flowing 
over  link  (i,k)  with  the  destination  j.  Note  that 
link  (k,i)  is  separate  from  link  (i,k)  and  it  is 
possible  to  have  > 0 and  > 0.  We 

assume  that  all  traffic  destined  for  j is  removed 
from  the  network  when  it  gets  to  j.  Thus  no  traf- 
fic for  j flows  outward  from  j and  is  de- 

fined only  for  i j,  (i,k)eJJ?  , tdiere  denotes 
the  set  of  links  in  the  network. 

We  assume  that  all  traffic  destined  for  j 
eventually  gets  there  by  travelling  over  paths  of 
the  network.  Thus  all  the  traffic  for  j that  comes 
into  each  node  i j must  go  out  of  it  again,  giv- 
ing us  the  fundamental  conservation  equation  for 
each  i ^ j , 


r^(j) 


£ 


(1) 


The  first  sum  edxive  is  over  the  integers  £ 
for  which  and  £ ><  j.  The  second  is  over 

the  integers  k for  which  (i,k)e5’;  for  brevity,  we 
regeurd  such  restrictions  as  understood  in  what 
follows.  Also  for  brevity,  we  shall  denote  the 
set  {rj^(j)}  of  inputs  over  all  i,j,  i j by  the 

vector  r (the  ordering  of  the  components  is  unim- 
portant) and  the  set  of  flows  for  all 

(i,k)sS^  and  all  j ^ i by  the  vector  f. 

Definition:  A set  of  flows  f ^ 0 is  a 
feasible  multicommodity  flow  for  the  in- 
puts r ^ 0 if  (1)  is  satisfied. 

We  regard  each  destination  j hero  as  cor- 
responding to  a commodity;  an  eunount  of  that 

commodity  enters  the  network  at  node  i and  is  moved 
to  j in  accordance  with  the  flow  vector  f.  Each 
such  feasible  multicommodity  flow  represents  a 
particular  way  of  routing  the  inputs  r to  their 
destination.  Such  routings  can  be  implemented,  in 
the  presence  of  dynamic  variations  on  the  inputs, 
by  taking  the  dynamically  varying  traffic  arriving 
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at  i dsstined  for  j emd  allocating  it  to  outgoing 
links  in  the  same  ratio  as  given  by  the  time  aver- 
age flows 

Sometimes  one  wishes  to  consider  ciulticoo- 
iDodity  flows  in  which  more  than  one  node  is  used 
as  a sink  for  a given  coinmodity.  This  situation 
can  be  reduced  to  the  one  just  described  by  add- 
ing an  extra  node  to  the  network  for  each  commodity 
and  adding  a link  from  each  sink  for  a commodity 
to  the  special  node  for  that  consnodity. 

For  purposes  of  studying  the  delay  experi- 
enced by  traffic  in  a network  (or  some  other  kind 
of  cost  experienced  by  the  cooimodities) , we  are 
often  Interested  primarily  in  the  aggregate,  or 
total  flows,  in  the  network  lines.  In  terms  of  a- 
multlconaodity  flow  f,  the  aggregate  flow  F with 
components  F^  is  given  by 

“ I ‘tvW  I (i.Wei?  (2) 

ax  ix 

Definition : F is  a feasible  aggregate 

flow  for  an  input  F if  (2)  is  satisfied 

for  some  feasible  multi-commodity  flow  f. 

There  appears  to  be  no  very  simple  test  for 
whether  or  not  F is  a feasible  aggregate  flow  and 
in  fact,  the  problem  is  almost  Identical  to  the 
following  classic  "multicommodity  flow  problem". 
Given  a set  of  capacities  ^ 0,  (i,k)e2’,  de- 
termine whether  there  is  a feasiole  multicomodlty 
flow  f for  which 


single  commodity  traffic  to  node  3 for  Figure  lb 
contains  a loop.  In  general,  the  uingle  commodity 
traffic  to  node  j contains  a loop  if  there  is  a 
cycle  of  nodes,  i,k,  1, . . . ,r,s,i  such  that 

0,  0,...,f^g(j)>  0,  f^^(j)  >0.  Loops 

of  this  type  have  been  treated  in  [1] . 


I » (i,k)e^ 


(3) 
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In  particular,  it  should  be  stressed  that  conserva- 
tion of  aggregate  flow,  as  expressed  by  (4)  below, 
is  necessary,  but  far  from  sufficient,  for  F to  be 

feasible. 
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-I 
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I r,(j)  - I r^(i), 
j m 


all  i 

(4) 


Loops  and  Loopfreedom  in  Network  Flows 


Figure  Ic  is  similar,  but  harder  to  handle, 
since  the  single  commodity  flow  to  destination  1 
is  loop  free  and  that  to  destination  2 is  loop 
free.  However,  the  same  aggregate  flow  is  achieved 
by  sending  the  traffic  from  3 to  2 directly  and 
that  from  3 to  1 via  node  4 with  a loop  to  2.  This 
leads  us  to  our  general  definition  of  loops  for 
multicommodity  flows. 


Before  defining  what  we  mean  by  loops  in  net- 
tK>rk  flows,  we  give  several  examples  in  Figure  1. 

To  avoid  cluttering  the  figures,  the  links  are  not 
shown  and  each  line  between  two  nodes  Indicates  a 
unit  flow.  The  flows  are  broken  down  into  individ- 
ual source-destination  flows  indicated  by  the  ex- 
tensions of  the  lines  through  the  nodes.  Figure 
la  shows  the  most  obvious  type  of  looping  where 
the  traffic  from  node  1 to  3 loops  from  1 to  2 and 
back  before  going  to  3.  Figiure  lb  shows  a situa- 
tion in  which  neither  the  source-destination  traf- 
fic from  node  1 to  3 nor  that  from  2 to  3 has  a 
loop.  However,  if  we  look  at  the  traffic  for  node 
3 as  a commodity,  we  have  ” ^23*^^  “ ^21*^^ 

“ ” !•  This  could  also  be  broken  into  indi- 

vidual source  destination  pairs  by  having  the  traf- 
fic from  source  2 to  3 go  directly  over  link  (2,3) 
and  the  traffic  from  source  1 to  3 go  over  the 
looping  path  of  Figure  la.  Thus  we  say  that  the 


Definition;  A feasible  aggregate  flow 
F for  the  input  r contains  loops  if  it 
can  be  expressed  as  F = F'  + F"  where 
F'  is  a feasible  aggregate  flow  for  r and 
F"  > 0 is  a conservative  flow.  F is  loop- 
free  if  it  does  not  contain  loops. 

A flow  F is  conservative  if  it  satisfies  (4)  with 
the  right  hand  side  set  to  zero.  Also,  by  F"  > 0 
we  mean  F"j^  ^ 0 for  all  (i,k)c2’  with  strict  un- 

etguality  for  at  least  one  (i,k)ei?  . He  shall  use 
the  terminology  positive  vector  to  refer  to  the 
stronger  condition  where  all  components  are  posi- 
tive. It  is  a well  known  result  of  graph  theory 
that  a conservative  flow  can  be  represented  as  a 
finite  cum  of  cycles,  or  loops.  ... 

Figure  Id  (due  to  J.  Hozencraft)  contains  no 
new  complications  and  is  included  toerely  to  Indicate 
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that  loo?rre<fdc.-3  ic  not  a prop-irty  that  car.  be 
checked  by  inspection.  Below  we  show  the  separa- 
tion of  P for  Figure  Id  into  F’  and  F".  Note  that 
simply  hawing  a loop  of  aggregate  flow  does  not 
destroy  loopfreedon;  consider  for  exxoiple  a 2 node 
network  with  traffic  flowing  both  ways. 
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Lemma  Is  A feasible  aggregate  flow  F contains 
loops  iff  F'  is  euiother  feasible  aggregate  flow 
for  the  same  inputs  with  F'  < F. 

Proof;  Suppose  F'  < F,  and  let  F"  = F - F*.  Then 
F"  < 0.  Since  F and  F'  each  satisfy  (4),  we  can 
subtract  these  equations  from  each  other,  getting 

I ^Lc-  ° 


called  a feasible  solution  to  the  L?p;  if  in  addi- 
tion they  minimize  (5)  they  are  called  an  optimal 
feasible  solution.  Note  that  if  f • 0 satisfies 

(7),  it  is  a feasible  multicommodity  flow  for  r. 
Thus  if  we  define  a vector  F'  with  components  F^ 

satisfying 

- I «^(j>  » a,k)e^,  (8) 

J t • 

we  see  that  F*  is  a feasible  aggregate  flow  for  r. 
Equation  (6)  can  then  be  rewritten 

s + y « F ■ (9) 

Note  that  if  OPT  <0,  then  there  is  a feasible 
solution  to  the  LPP  with  s >0,  and  from  (9),  F 
contains  loops.  Conversely,  if  OPT  - 0,  there  is 
no  s > 0 and  feasible  F'  satisfying  (9)  and  F is 
loopfree.  Next  note  that  the  LPP  always  has  at 
least  one  feasible  solution,  obtained  by  setting 
s s 0 and  setting  f equal  to  the  multlcocmodity 
flow  corresponding  to  F.  Since  s is  bounded, 

0 s < F,  the  LPP  must  have  an  optimal  feasible 
solution  (although  generally  not  unique) . Finally, 
from  (9)  and  (5),  we  see  that  the  F*  in  any  optimal 
feasible  solution  is  loop  free.  Summarizing,  we 
have  the  following  theorem. 

Theorem  1;  Let  F be  a feasible  aggregate  flow  for 
the  input  r.  Then  the  minimum  in  (5)  exists  and 
the  corresponding  F*  is  loopfree;  F is  loopfree  iff 
OPT  = 0. 


Thus  F"  is  conservative  and  F contains  loops.  Con- 
versely, if  F contains  loops,  there  is  an  F’  < F 
by  definition. 

A Linear  Progr^Ulmling  Approach  to  Loopfreedom 


In  this  section  we  formulate  a linear  program- 
niing  problem  to  determine  whether  or  not  a feasible 
aggregate  flow  is  loopfree.  Our  purpose  is  not  to 
develop  a computational  procedure  but  rather  to  use 
the  known  results  of  linear  programming  to  deter- 
mine some  of  the  consequences  of  loopfreedom.  Let 
r be  the  input  to  a given  network  and  let  F be  a 
given  feasible  aggregate  flow.  Let  s be  a vector 
of  L slack  components  with  one  component  s^j^  for 

each  (i,k)e2’.  Let  f*  be  a vector  of  (n-l)L  com- 
ponents with  one  component  fot’  each  (i,k)c2? 

j i.  The  linear  programming  problem  (LPP)  is  to 
choose  s ^ 0,  f ' ^ 0 to  fora  the  minimization; 

OPT  ” min  J “ {5) 

(i,k)eS> 


subject  to  the  following  L + n(n-l)  linearly  inde- 
pendent constraints 

‘ik  j " ^ik  ' 

I f^(j)  - I £;,(j)  “ r,<j)  » all  i,j,  i / j 

ik  /.i  i 

(7) 


Next  we  derive  the  dual  linear  programming 
problem  to  that  in  (5)  to  (7) . The  standard  form 
for  a LPP  is  to  choose  x ^ 0 to  minimize  c’x  sub- 
ject to  A*x  = b for  a given  row  vector  c,  matrix  A, 
and  column  vector  b.  The  dual  LPP  to  this  stzmdard 
primal  LPP  (see,  for  example,  Luenberger  (21)  is  to 
choose  a row  vector  u to  maximize  u'b  subject  to 
u’A  ^ c.  The  duality  theorem  states  that  if  both 
the  primal  and  dual  have  feasible  solutions,  or  if 
either  has  an  optimal  feasible  solution,  then  both 
have  optimal  feasible  solutions  with  the  same  value. 
In  our  case,  since  the  primal  always  has  an  optimal 
feasible  solution,  the  dual  does  also. 

In  order  to  visualize  (S)  to  (7)  in  standard 
fora,  let  X be  made  up  of  the  L components  of  s fol- 
lowed by  the  (n-l)L  components  of  f,  thus  x has  nL 
components.  Choosing  c as  L components  of  value  -1 
followed  by  (n-l)L  zeros,  (5)  is  equivalent  to  min- 
imizing c'x.  It  turns  out  that  the  dual  variables 
will  have  more  significance  If  we  multiply  (6)  on 
both  sides  by  -1.  Then  the  vector  b will  be  the 
components  of  -F  followed  by  r. 

Finally  we  must  associate  a dual  variable  with 
each  equation  in  (6) , (7) . Let  d^^j^  be  associated 

with  (6)  for  each  (i,k)ci?,  let  D^(j)  be  associated 

with  (7)  for  each  i J,  and  let  d,  D be  vectors 
with  these  components  respectively.  Letting  u bo 
the  components  of  d followed  by  those  of  D,  the  dual 
problem  is  to  form  the  maximum; 


If  s ^ 0, 


f ^ 0 satisfy  (6)  and  (7), 


they  are 
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OPT 


ISOX 

d,D 


<i,k)e5?  i.3  ^ ^ 

If'j 


(10) 


a miniaua  distance  path  fron  i to  j.  Similarly  vi» 
say  that  ? is  a shortest  route  aggregate  for 

d if  some  Dulticomaodity  flow  f corresponding  to  F 
is  a shortest  route  multicommodity  flow  for  d. 


subject  to  uA  ^ c.  This  constraint  can  be  viewed 
as  nL  inequalities > one  for  each  column  of  A,  or 
equivalently > one  for  each  component  of  x.  Iden- 
tifying the  elements  of  A with  the  coefficients  of 
a and  f in  (6)  and  (7),  these  nL  inequalities  can 
be  written  explicitly  for  each  (i,k)ti2?  as 


“ih  - " 


0,(j)  1 t D^(j) 


each  j f>  i,k 


(11) 

(12a) 


We  now  state  the  complementary  slackness 
theorem  of  linear  progranning  in  the  context  of  the 
primal  and  dual  problems  here  and  then  show  its 
relation  to  shortest  route  flows. 

Theorem  2;  Let  s,  f be  a feasible  solution  to  the 
primal  problem  and  d,  D be  a feasible  solution  to 
the  dual.  A necessary  and  sufficient  set  of  con- 
ditions for  both  solutions  to  be  optimal  is  for 
all  (i,k)e^,  i j. 


I 


L 


I 


Dj,(k)  < d^  (12b) 

Inequality  (10)  corresponds  to  the  variable  s^j^, 

(12a)  to  the  variables  fjj^(j)»  a“d  (12b)  to 

To  avoid  the  notational  inconvenience  of  constantly 
sep^trating  the  case  j * k from  j ^ k in  (12),  we 
shall  often  use  (12a)  for  both  cases  under  the  con- 
vention that  Dj ( j)  is  a constant  equal  to  zero. 

In  order  to  understand  the  meaning  of  (12) , 
suppose  that  a path  of  links  from  node  i to  j con- 
sists of  the  links  (i,k) , (k,!') , . . . , (m,  j) . By  .ap- 
plying (12)  recursively  to  the  final  element  in 
(12a),  we  get 


°i‘l>  1 ^ik^-Scl*  •••  ■^‘'mj 

If  we  interpret  d^  as  a length  associated  with 

link  (l,k),  and  take  the  distance  on  a path  to  be 
the  sum  of  the  link  lengths  on  the  path,  then  (13) 
states  that  D. (J)  is  less  than  or  equal  to  the 
distance  on  each  path  from  i to  j for  the  given  d. 
If,  for  a given  d,  we  simply  choose  each  D^(j)  as 

the  distance  of  the  minimum  distance  path  from  i 
to  j (according  to  the  lengths  d) , then  this  D 
clearly  satisfies  (12) . From  (13) , any  other 
choice  D*  will  satisfy  D*  < D.  The  quantity  (10) 
to  be  suucimized  in  the  dual  problem  multiplies  each 
D^(j)  by  the  non-negative  coefficient  r^(j),  estab- 
lishing the  following  leissa. 


Lemma  2:  For  given  positive  d,  (10)  is  maximized 
over  D,  subject  to  (12)  by  setting  D^(j),  for  each 

l,j,  i^j,  equal  to  the  minimum  distance  (i.e., 
the  length  of  the  minimum  distance  path)  from  1 to 
j for  the  given  d. 

Minimum  distance  paths  are  closely  related  to 
the  concept  of  shortest  route  flows.  Specifically, 
for  a given  set  of  positive  link  lengths  d and  the 
corresponding  minimum  distances  D we  say  that  a 
multiconmodity  flow  f is  a shortest  route  nulti- 
coCTnodlty  flow  for  d if  ^ whenever 

D^(j)  < d^jj(J);  i.e.,  whenever  link  (i,k)  is  not  on 


Sjjj  - 0 if  d^^  > 1 / (i,k)e^  (14) 

f|j^(j)  - 0 if  Dj^(j)  < d^^  + D^(j)  (15) 

The  following  theorem  relates  loop  freedom  to 
shortest  route  flows. 

Theorem  3;  Let  F be  a feasible  aggregate  flow.  F 
is  loopfree  iff  there  is  a set  of  positive  link 
lengths  d such  that  F is  a shortest  route  aggregate 
flow  for  d.  Furthermore,  any  such  d,  scaled  up  so 
that  its  components  exceed  one,  yields  an  optimal 
solution  to  the  dual  problem,  (10)  to  (12) . 

Proof;  Let  f be  any  multicommodity  flow  correspond- 
ing to  F,  and  note  that  s « 0,  f « f yields  a 
feasible  solution  to  the  primal  problem.  Assume 
first  that  F is  loop  free,  so  that  the  alxive  solu- 
tion is  optimal.  Choose  link  lengths  ^ 1 and 

the  associated  minimum  distances  D to  give  etn  op- 
timal solution  to  the  dual.  From  Theorem  3,  (15) 
is  satisfied  and  hence  f = f is  a shortest  route 
flow.  Next  assume  that  f is  a shortest  route  flow 
for  some  set  of  positive  link  lengths  d* . Let 
d,.  « d!  /min  dl  for  all  (i,k)e^.  Thus 

d^j^  ^ 1 and  f is  a shortest  route  flow  for  d. 

Since  d,  D is  a feasible* solution  to  the  dual,  and 
since  (14)  and  (15)  are  satisfied,  we  have  an  op- 
timal feasible  solution  to  the  primal  and  F is 
loopfree. 

Part  of  the  reason  for  our  Interest  in  this 
result  stems  from  our  conjecture  that  in  good 
quasi-static  routing  algorithms,  the  act  of  chang- 
ing a flow  f^jj(j)  from  zero  to  non-zero  or  vice 

versa  (i.e.  establishing  or  terminating  routes) 
should  be  done  much  more  slowly  and  carefully  than 
the  act  of  changing  a flow  between  different  pos- 
itive values.  Since  the  question  of  whether  or 
not  a multicommodity  flow  is  a shortest  route  flow 
(for  a given  d)  is  simply  a question  of  which  flows 
are  zero,  it  appears  that  these  lengths 
should  play  a role  in  quasi-static  algorithms. 

Our  next  few  results  concern  the  basic  optimal 
solutions  of  the  primal  and  dual  LPP’s.  The  matrix 
A in  these  problems  is  an  L t n(n-l)  by  Ln  matrix 
with  linearly  independent  rows.  A matrix  B is 
called  a basis  for  the  LPP  if  it  consists  of 


5 


L + n(a-l)  linearly  independent  columns  or  A (i.e., 
it  is  formed  from  A by  deleting  linearly  dependent 
columns) . An  x satisfying  Ax  = b is  called  a basic 
solution  (with  basis  B)  if  the  components  of  x 
corresponding  to  columns  of  A not  in  B are  all 
zero.  If  we  let  x be  the  vector  of  components  of 
X corresponding  to® the  columns  in  3,  then  we  see 
that  Xg  is  uniquely  determined  by  Bx^  “ bs  thus 

there  is  one  basic  solution  x corresponding  to  each 
basis  B.  If  this  basic  x also  satisfies  x > 0,  it 
is  called  a basic  feasible  solution.  Another  well 
known  result  of  linear  prograimning  (21  (and  the 
justification  for  the  simplex  algorithms)  is  that 
if  (optimal)  feasible  solutions  exist  to  the  primal 
then  a Isasic  (optimal)  feasible  solution  exists. 

Basic  dual  solutions  are  defined  in  a similar 
way.  Let  c be  the  vector  of  components  of  c cor- 
responding ^ the  columns  of  A in  B.  Then  u is  a 
basic  solution  to  the  dual  problem  if  it  satisfies 

uB  a c . Again  u is  unique,  given  B,  and  u is  a 
B 

basic  feasible  solution  if  it  also  satisfies 
uA  ^ c. 

Our  problem  now  is  to  investigate  what  sets 
of  columns  of  A form  a basis  B.  The  simplest  ap- 
proach is  in  terms  of  the  dual  constraints  (11)  and 
(12).  We  must  find  a subset  of  L + n(n-l)  of  these 
inequalities,  idiich  when  satisfied  with  equality, 
uniquely  specify  d,  D;  such  a subset  is  called  a 
set  of  basic  equations,  and  of  course  maJtes  up  a 
basis  B.  First  note  that  for  each  destination 
commodity  j,  (12)  consists  of  a set  of  inequalities, 
one  for  each  lin)c  not  originating  at  node  j.  These 
are  the  only  inequalities  that  involve  the  n-1  vari- 
ables Dj^(j),  i j,  and  thus  at  least  n-1  of  them 

must  be  used  as  basic  equations.  Next  observe  that 
the  basic  equations  from  (12)  for  a given  j must 
correspond  to  a set  of  links  T^  that  contain  a 

spanning  tree  of  the  network  (where  the  orientation 
of  the  links  in  the  tree  is  immaterial) . Other- 
wise there  %e3uld  be  a nonempty  set  of  nodes  that 
were  unconnected  to  j by  the  links  of  T^,  and  the 

equations  corresponding  to  T^  could  not  uniquely 
specify  D^(j)  for  the  i in  the  set.  The  equations 

corresponding  to  this  spanning  tree  for  a given 
destination  J can  now  be  rewritten  to  express  each 
D.  (j)  as  a sum  of  signed  link  distances  from  i to 
j , using  the  path  in  the  tree  from  i to  j . If  the 
link  (k,!)  on  the  path  from  i to  j has  the  same 
direction  as  the  path,  then  is  used  in  the 

sum  and  otherwise  is  used. 

He  have  now  identified  n(n-l)  of  the  basic 
equations  that  must  be  used  to  uniquely  specify  d, 

D (i.e.  for  each  of  n destinations,  there  are  n-1 
equations  corresponding  to  the  links  of  a spanning 
tree) . These  equations  uniquely  specify  D as  a 
function  of  d.  Some  destinations  might  have  more 
than  n-1  of  their  inequalities  used  as  basic  equa- 
tions. Each  additional  such  equation  D^(j)  • 

♦ Dj^(j)  can  be  written  using  our  previous  solu- 
tion for  D^(j)  and  D|j(i)>  this  equation  now  becomes 


tlie  sum  of  signad  link  distances  around  a cycle 
set  equal  to  zero.  This  is  the  same  type  of  equa- 
tion as  one  uses  in  an  electrical  network,  setting 
the  sum  of  the  voltages  around  a cycle  equal  to 
zero.  As  is  well  Icnown  from  electrical  circuit 
theory  (and  from  graph  theory) , there  exist  at  most 
L - ntl  linearly  independent  such  equations,  and 
since  these  equations  Involve  only  d,  this  is  a 
maximum  over  all  destinations  together.  We  con- 
clude from  this  that  at  least  n-1  of  the  inequal- 
ities d^  _>  1 must  be  included  as  equalities  in 

any  basis.  Although  we  do  not  need  it  here,  the 
links  corresponding  to  these  equalities  must  also 
include  a spanning  tree. 

Finally,  suppose  we  choose  a basis  with  n(n-l) 
spanning  tree  equations  and  exactly  L - n+1  inde- 
pendent cycle  equations.  Since  these  cycle  equa- 
tions are  homogeneous  and  form  a basis  for  the  set 
of  cycle  equations,  the  sum  of  the  distances  around 
every  cycle  must  be  zero;  thus  the  dual  solution 
from  any  such  basis  is  not  feasible.  We  summarize 
these  results  in  the  following  theorem. 

Theorem  4t  Every  basis  for  the  LPP  of  (5)  to  (7) 
contains  at  least  n-1  columns  corresponding  to 
slack  variables  s^j^  and,  for  each  j,  there  is  a set 

T.  of  at  least  n-1  linics,  which  include  a spanning 

tree,  such  that  the  column  corresponding  to 

for  each  (i,k)  eT^  is  in  the  )3asls.  xf  the  basis 

is  feasible  for  the  dual,  there  are  at  least  n 
columns  corresponding  to  slack  variables. 

Next  we  consider  the  restrictions  on  a basis 
B if  the  corresponding  primal  basic  solution  is  to 
Ise  feasible.  First  assume  that  r^(j)  >0  for  all 

i,j,  i j.  Then  every  feasible  solution  to  the 
primal  isust  contain,  for  each  j,  a directed  span- 
ning tree  to  j of  links  for  which  f]jj(j)  ^ 0.  By 

a directed  spanning  tree  to  j we  mean  a tree  such 
that  each  node  i has  a directed  path  to  j through 
the  tree  links;  i.e.  a tree  in  which  all  links 
point  toward  j.  Since  an  optimal  basic  feasible 
solution  to  the  primal  must  )iave  a basis  for  which 
the  dual  is  also  feasible,  we  see  that  the  basic 
variables  for  such  a basis  include  at  least  n slack 
variables,  and  for  each  j,  a directed  tree  to  j of 
variables  Because  of  the  directed  trees  in 

the  basis,  it  is  easy  to  see  that  the  corresponding 
optimal  basic  dual  solution  d,  D has  the  property 
that  the  components  of  0 are  minimum  distances  cor- 
responding to  d. 

Theorem  5t  Let  F be  a feasible  aggregate  flow  for 
input  r ^ 0.  In  any  optimum  feasible  solution  to 
(5)  to  (7) , F*  (and  F if  F is  loopfree)  is  a short- 
est route  flow  for  some  d such  that  (d,D)  is  a 
basic  feasible  solution  to  the  dual  2md  the  com- 
ponents of  D are  minimum  distances  for  d. 

Proof;  We  have  just  seen  that  if  r^^d)  > 0 for  all 

i>j,  i 7^  j>  then  there  is  a iMsic  optimal  feasible 
solution  (d,D)  to  the  dual  for  which  the  components 
of  D are  minimum  distances  for  d.  From  Theorem  2, 


F‘  (and  r if  F is  loopfrae)  is  a sliortos:  route 
flow  for  d in  any  optiinal  solution  of  the  prinal. 

If  r^d)  » 0 for  one  or  more  i,j,  i j,  we  let 

r(C)  be  an  input  formed  by  adding  e>  0 to  each 
r^(j)  = 0.  Similarly  let  when 

for  each  (l,k)  e^is  chosen  to  make  F(E)  a fea- 
sible aggregate  flow  for  r(c).  Let  c vary  over  a 
sequence  -*•  0.  Since  there  are  a finite  number 

of  basic  feasible  solutions  to  the  dual,  and  these 
are  independent  of  r,F,  we  see  that  one  of  these 
oust  satisfy  Theorem  5 for  am  infinite  set  of  e^, 

auid  thus  must  satisfy  Theorem  5 in  the  limit  C=>  0. 

Theorem  4 and  5 have  an  important  consequence 
with  respect  to  the  amount  of  alternate  routing 
required  to  achieve  a given  F (or  its  looafree  re- 
duction F*).  Let  0^(j)  be  the  number  of  links 

(i,k)  for  which  > 0 for  a given  basic  opti- 

mal solution  to  the  primal,  amd  let  a be  the  aver- 
age of  over  the  n(n-l)  node  pairs  i,j,  i ^ j. 

Then  a • 1 corresponds  to  no  alternate  routing. 

He  have  just  shown  that  in  addition  to  the  directed 
spanning  trees  (which  contribute  1 to  each  a^(j)), 

basic  optimal  solutions  have  at  most  L - n addi- 
tional non  zero  over  all  i,k,j.  This  es- 

tablishes the  following  theorem. 

Theorem  6;  The  amount  of  alternate  routing  in  a 
basic  optimal  solution  to  the  I.FP  of  (5)  to  (7)  is 
bounded  by 


Since  L is  typically  much  smaller  than  n(n-l) 
in  luge  net%«orks,  very  little  alternate  routing  is 
required  to  achieve  any  desired  loopfree  aggregate 
flow. 

In  order  to  obtain  a better  understanding  of 
the  meaning  of  Theorem  5,  let  be  the  sec  of  all 

feasible  aggregate  flows  F for  a given  input  r, 
and  let^^(d)  be  the  set  of  shortest  route  aggre- 
gate flows  for  input  r and  link  length  assignment 
d.  Let^  be  the  set  of  link  length  assignments  d 
for  which  (d,D)  is  a basic  feasible  solution  to 
(10)  to  (12)  and  for  which  D is  the  set  of  minirum 
distances  for  d.  Theorem  5 then  asserts  that  the 
set  of  loopfree  flows  for  r,  LFF^,  as  given  by 

LFF  - 9 (d)  (17) 

' dtss 

He  can  Interpret  .:^(d)  (for  any  d with  posi- 
tive components)  as  the  set  of  vectors  F that  min- 
imtimlze  F-d  over  the  constraint  that  F be  a feasi- 
ble aggregate  flow  for  r (this  comes  directly  from 
the  definition  of  a shortest  route  aggregate  flow) . 
He  also  obtain  the  geometric  interpretation  from 
Figure  2 that  if  F minimizes  F*d  over  F e JF^, 

t)ian  any  vector  F'  for  which  F*  < F must  satisfy 
F'*d  < F*di  thus  F'  < F implies  F'  i S =o  that 


riOM) 


^ • (ri  r < r«) 


Caeaatrle  tatwpMtAtloii 
rifur*  3 


F is  loopfree  (a  fact  we  already  )cnow  from  Theorem 
3).  More  compactly,  this  says  that  ^^(d)  is  a 

face  of  SF  , i.e.  the  subset  of  ^ that  lies  on 
r r 

the  supporting  hyperplane  of  y with  normal  vector 

d. 

Theorem  7;  Let  d be  an  arbitrary  positive  length 
vector.  Then  ,y^(d)  C ^^(d’)  for  some  d*  c ^ • 

He  postpone  the  proof  of  this  until  after  the 
proof  of  Theorem  8.  The  theorem  says  that  all  the 
shortest  route  flows  for  any  given  positive  d are 
included  in  the  shortest  route  flows  for  one  of  the 
basic  distance  vectors  in  ^ . Geometrically,  it 
says  that  the  faces  ^^(d)  for  d are  in  some 

sense  superfluous,  since  each  of  them  is  contained 
in  some  basic  face  ^ (d').  It  should  also  be 
noted  from  the  dual  problem  (10)  to  (12)  that  the 
set  of  basic  distance  vectors  is  Independent  of 
the  input  r.  . Thus  as  r changes,  each  of  the  basic 
faces  y (d)  move  in  space  but  maintain  the  same 


Some  additional  interpretation  of  these  faces 
arises  from  defining 


{(i,k,j)t  [r^(j) 


dik^Dkd)} 


where  D is  the  set  of  minimum  distances  for  d. 
^^(d),  then,  is  the  set  of  aggregate  flows  for 

which  a corresponding  multicommodity  flow  exists  in 
which  each  commodity  j flows  only  over  links  (i,k) 
for  which  (i,k,j)  is  in  .‘^(d).  This  means  that 
^^(d)  is  a function  of  d only  through  the  set 

^ (d) . It  turns  out  that  several  basic  distance 
vectors  can  have  the  same  ^(d)  (in  which  case  they 
determine  the  same  face  JF  (d) ) . It  can  also  hap- 
pen that  ^(d)  for  one  basic  d can  be  strictly  con- 
tained in  J?(d’)  for  another  basic  d',  in  which  case 


is  strictly  contained  in^(d')  (assicaing 

positive  r) . Thus,  in  some  sense,  a more  elegent 
representation  of  the  class  of  loop  free  flows 
could  be  obtained  directly  in  terms  of  a minimal 
class  of  sets  SUd) , but  we  have  been  unable  to  do 
this  except  through  the  basic  distance  vectors. 

Theorem  8:  Let  F. . F be  feasible  aggre- 

“ 12  n 

gate  flows  for  input  r and  let  A,,  X,,...,  X be 

X z m 

positive  numbers  summing  to  1.  Then  “ 

j X^  Fj^  is  in^(d)  iff  F^  e ^^(d)  for 

i-i  - - 

1 < i < m. 

Discussion;  This  says  that  the  convex  combination 
of  loop  free  flows  is  not  in  general  loopfree  un- 
less each  of  the  combined  flows  are  shortest  route 
flows  for  the  same  d.  This  theorem  appears  to  give 
an  indication  why  routing  algorithms  such  as  the 
Cantor-Gerla  algorithm  [31 , which  operate  by  tak- 
ing convex  combinations  of  extremal  flows  (i.e. 
flows  that  2tre  extremal  points  of  the  set  , con- 
verge quickly  at  first  and  then  more  slowly.  Each 
new  extremal  flow,  when  combined  with  the  others, 
typically  adds  a new  set  of  loops  to  the  solution, 
and  these  loops  are  eliminated  only  when  all  of  the 
extremal  flows  being  used  are  shortest  route  for 
the  same  d. 

Proof  of  Theorem  8;  Let  c “ min  F*d  over  F . 
n 

Fn’d  - I X F .d  > c 
i-'l 

with  equality  iff  =■  c for  1 i £ m. 

Proof  of  Theorem  7 ; Suppose,  to  the  contrary  that 
is  not  contained  in,y^(d')  for  any  d' 

Let  d, , d.,...,  d be  the  elements  of  ^ and  let  F. 

1 2 IQ  X 

be  an  aggregate  flow  in  ^ (d)  but  not  in  ^ (d^) ; 

the  F,  need  not  be  distinct.  Let  X.,...,  X be 
1 1 m 

positive  numbers  summing  to  1,  and  let  F • 


References 

R.  Gallager,  "A  Minimum  Delay  Rout- 
ing Algorithm  Using  Distributed  Com- 
putation", IEEE  Trans.  Commun. , Vol. 
COM-25,  No.  1,  Jan.  1977,  pp.  73-85. 

D,— G»— tuenberger^,  - lHfero4uct-i<«— to 
Linear  and  Non-Linear  Programming, 
Reading,  Ma.,  Addison  Wesley,  1973. 


D.  G.  Cantor  and  M.  Gerla,  "Optimal 
Routing  in  a Packet  Switched  Computer 
Network",  IEEE  Trans.  Comput. , Vo.  j 
C-23,  pp.  1062-1069,  Oct.  1974. 


J From  Theorem  8,  F.  is  not  in  ^ Id.), 

i_j  1 1 o r 1 

for  1 £ 1 £ m,  and  from  (17)  is  not  a loopfree 
flow.  However,  again  from  Theorem  8,  F^  is  in 
^(d),  and  hence  is  a loopfree  flow,  supplying  the 
desired  contradiction. 
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